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(54) 3D seismic surveying 

(57) A method for determining a representation of an anisotropic earth formation is described, wherein the 
Inversion step for generating the representation from recorded seismic data is performed using several values 
of a pre-defined anisotropy parameter addition to components of the moveout velocity. The anisotropy 
parameter is preferably introduced into the moveout equation or dispersion relation, which in turn is used in 
the inversion process. A more Accurate representation of earth formation can be generated. 
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Inversion Method for Seismic Data 

The present invention relates to an inversion or migration 
method for seismic data. More specifically it relates to such a 
5 method as applied to seismic data characterizing anisotropic 
media. 

BACKGROUND AND OBJECT OF THE INVENTION 

10 

It is well recognized that most rock formations in the Earth are 
to some extent anisotropic. For the most part, the anisotropy 
has been ignored in seismic processing. Its effects are small 
when compared to assumptions such as plane-layer approximation 

15 imposed by NMO and stack processing or the 2D earth model 
approximation imposed by 2D processing. However, with recent 
advances in inverse problem theory and supercomputing 
technology, most of these approximations have been overcome, 
leading to the development of 3D prestack imaging and inversion. 

20 Recent studies suggest that the type' of details and accuracy 
expected from 3D prestack imaging and inversion can be distorted 
by neglecting anisotropy. 

To include anisotropy in 3D prestack imaging and inversion, we 
25 need an accurate estimation of the anisotropic background 
velocity. The reconstruction of anisotropic background velocity 
models is probably the most important challenge in processing 
with anisotropy. Conventional methods for reconstructing the 
background velocity model, like migration-velocity methods, 
30 often assume an isotropic subsurface and can yield to inaccurate 
reservoir descriptions when the subsurface contains anisotropic 
rock formations. 



The knovm migration-velocity method consists of scanning over 
different velocity models using prestack migration, the "good" 
velocity model is then constructed based on focusing analysis. 
Further details of this method are described for example by O. 
Yilmaz, Seismic Data Processing, Society of Exploration 
Geophysists , 1987 . 

An inversion method for data in the space-time (x-t) domain has 
for example been described by A. Tarantola in: Geophysical 
Prospecting. No. 32 (1984), 998-1015). 

Another inversion process for data in the wavenumber- frequency 
(co-k) domain has for example been described by L. Ikelle et al . 
in: Geophysics, Vol. 51, No. 6 (June 1986), 1266-1276. 

An inversion method for isotropic media, also including a 
background velocity reconstruction, is described by L. Ikelle 
in: Geophys.J.Int . , No. 123(1995), 507-528. 

An attempt to reconstruct a velocity model for a transversely 
isotropic medium is described by T.A. Alkhalifah and I. Tsvankin 
in: 64th Ann. Intern. Mtg., Soc. Expl. Geophys . (1994), Expanded 
Abstracts, 1000-1003 . 

One reason why the isotropic linearized inversion or migration 
in the (co-k) domain are widely used through the industry for 
reconstructing background velocity models is that the dispersion 
relationship which controls the extrapolation operation is 
explicit. This aspect has been extensively exploited over the 
last two decades to produce efficient migration and linearized 
inversion algorithms in the (ci>-k) domain. For anisotropic media, 
the dispersion relationship is no longer explicit, it requires 
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numerical resolution of resolution of an eigenvalues-eigenvector 
system for each point in the (co-k) domain. Such computation is 
very slow. It renders inversion and migration, therefore 
anisotropic background reconstruction, economically 

5 unattractive. 

In view of the above cited prior art it is an object of the 
invention to improve the known inversion methods for seismic 
data such as to include anisotropy in the subsurface formation. 

10 

SUMMARY OF THE INVENTION 

The above-mentioned objects are achieved by methods as set forth 
15 in the appended claims. 

An important feature of the invention is that the inversion step 
for generating the representation from recorded seismic data is 
performed using several values of a pre-defined anisotropy 
20 parameter addition to components of the moveout velocity. The 
anisotropy parameter is preferably introduced into a moveout 
equation or dispersion relation, which in turn is used in the 
inversion process. 

25 The term inversion is defined to migration methods for seismic 
data, in particular any form of velocity analysis of such data. 
The data are assumed to cover the surveyed foinnation in more 
than one direction on the surface (3-D survey), preferably in 
the form of multi-offset and multi-azimuth data. 

30 

A preferred pre-conditioning step includes the removal of direct 
waves and multiple reflections from the data. However a method 



in accordance with the invention is best applied to prestack 
data, i.e., before the single recorded traces are added or 
otherwise combined. 

The invention can be described as generalizing the migration- 
velocity concept in two aspects: (i) the migration is replaced 
by multi-parameter linearized inversion for common azimuthal 
sections and (ii) the background velocity can be anisotropic in 
its elastic behavior. The combination of both improvements leads 
to a new scanning procedure for reconstructing anisotropic 
background media. 

It is regarded as a further important feature of the anisotropic 
or anelliptic parameter, which are preferably introduced into 
the dispersion or moveout equations, are dimensionless . 
Dimensionless is defined for the scope of the present invention 
as having the nature of a pure number, with a value independent 
from the choice of units. 

These and other features of the invention, preferred embodiments 
and variants thereof, and advantages will become appreciated and 
understood be those skilled in the art from the detailed 
description and drawings following hereinbelow. 



BRIEF DESCRIPTION OF DRAWINGS 



FIGS. lA, IB show a dispersion relationship in accordance with 
an example of the invention in two directions. 



FIG. 2 shows an anisotropic earth model- 
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MODE(S) FOR CARRYING OUT THE I^A^NTION 

The detailed description of an example in accordance with the 
5 invention, following hereinbelow, refers firstly to a method of 
deriving an anisotropic dispersion relationship which can be 
used in a linearized inversion or migration process. Secondly, 
an example of reconstructing the background or strata velocity 
model is described based on the anisotropic dispersion relation. 

10 

Although most rocks in the Earth are anisotropic, fortunately 
they are weakly anisotropic. Explicit formulae of the dispersion 
relationship even in anisotropic media can thus be derived by 
using the perturbation techniques with the isotropic component 

15 of the medium as the unperturbated mediLim. The perturbation 

method consists of expanding the eigenvalue-eigenvector system 
into perturbation series and keeping only the leading terms, 
• preferably the first two. The first term corresponds to the 
isotropic solution which is well Icnown. The second terms 

20 represents a first order perturbation from the isotropic 

solution. In the following an explicit representation of this 
second term is presented. 

If k^, kj are the horizontal wavenumbers corresponding to x and y 
25 directions, respectively, and if CO is the angular frequency 

corresponding to time t, the relationship between the vertical 
wavenumber, denoted q(z) and the angular frequency CO is called 
the dispersion relation. In the isotropic case for example, the 
dispersion relationship for P-wave is 



30 
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Herein QjCz) denotes the vertical wavenumber in the isotropic 
media; q(z) is the reserved symbol for the vertical wavenumber 
5 in the weakly anisotropic media. This convention is followed 
through this description. The vertical wavenumber q.,(-z) is 
dependent of (D, k^, k, and z; this holds also for q(z) . 

The vertical wavenumber q(z) is a solution of an eigenvalue- 
10 problem for which an analytic solution, q,(z), is known if the 
medium is isotropic. In the weak anisotropic case, q(z) deviates 
slightly from its underlying isotropic counterpart q^i^) - So, by 
using the perturbation technique, we can determine the vertical 
wavenumber q{z) in the form 



15 



25 



[2] q(z) = ^[q^Cz) f + 2qo(z)6q(z) 



where 5q(z) is the first order correction to cijz) . As there is 
no unique way of determining 5q(z) , the present example seeks a 
20 solution which provide the right trade-off between the level of 
accuracy and the algorithmic efficiency. The asymptotic 
expansion is used to develop q{z) as an asymptotic series in the 
density normalized stiffness tensor, a{ijpq). Thus, 6q(z) can be 
obtained as a linear combination of Aa{ijpq}. 



Another important form of the dispersion relationship consists 
of expressing the frequency 0) in terms of q, and k,. It is 
useful in migration and inversion with constant velocity. The 
frequency G) is decomposed into 
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[31 



CO = cOq + 6co 



■ where (O^ is the frequency associated with the oscillations in the 
5 isotropic component of the weakly anisotropic media and 8q) is the 



first order correction to co, due to the anisotropic component of 
the weakly anisotropic medium. The frequency (a, can be determined 
in terms of q,, k, and k, directly from equation [1] . To determine 
8©, the phase slowness has to be corrected as well. Let the 
10 unperturbed phase slowness be Yo and the actual phase slowness be 
y _ Yo+5y at a given wavenximber q,. Then 



15 which leads to the following relationship between the slowness 
and frequency perturbations at a given wavenumber 



[4] 



(5) 




20 or equivalently to 



[6] 




So far 5q can be the perturbation of the vertical wavenumber of 
25 qS-waves or qP-waves. For the remainder of this description, the 
example focuses on qP-waves only. However, an analysis similar 



- 8 - 

to that described below for qP-waves can be carried out for qS- 
waves by using the qS-wave expression of 5q. 



10 



15 



It should be also noted that the approximation (6) to the 
dispersion relationship is different from those proposed by 
Dellinger et al . in: J. of Seismic Explor. , No. 2(1993), 23-40. 
Contrary to those approximations, the approximation of eq. [6] 
is not restricted to transversely isotropic (TI) media and more 
importantly it can be used for prestack F-K migration. 

The vertical wavenumber q(2) for a transversely isotropic media 
can be written in the following form 



t7,q(z) - --^1 - 3 ^4 

V 



When T|=0 equation [7] can be written in elliptical form: 



[8] 



20 Therefore T| is the anelliptic or anisotropy parameter in the 

equation [7]. The ellipticity of this equation is described by 
the normal moveout velocity and the vertical velocity V^. The 
anisotropy parameter, the normal moveout velocity and the 
vertical velocity are important parameters in the scanning 

25 procedure for reconstructing anisotropic background velocity 

models as is described below. As for the following description, 
the term 'anelliptic* and "anisotropic* both denote directional 
dependency of physical properties of an earth formation. Thus, a 
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relation between the anisotropy or anelliptic parameters and the 
elastic coefficients of the stiffness tensor exists. 

In the following the dispersion relationship, through the 
5 vertical wavenumber q(2), for an azimuthally anisotropic media 
is derived. In the present example, an orthorhombic medium is 
considered. The orthorhombic case is a standard type of 
anisotropy and it ,is often encountered in the Earth exploration. 
The dispersion relation is 

10 



q(z) = — X 
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15 If the anisotropy parameters (i=4,5,6) all equal 0, equation 
[9] describes an ellipsoid. It becomes: 



[10] 



.NMO 1 

CO 



+ 



V k^ 

V.NMO 2 

CO 



= 1 



20 So, the anellipticity in the orthorhombic media is described by 
the three parameters . The ellipsoidal behavior is described by 
a vertical velocity and two elliptical normal moveout 
velocities V, and V^. These parameters will be used in the 
scanning procedure described below. 



It is important to notice that the above definition of the 
anisotropy parameters is not unique. An alternative definition 
is for example given by multiplying each parameter as defined in 
eq. [9] by factors of f(V^/V,), where f is an arbitrary function 
of the argument. This redefinition or scaling would effectively 
eliminate the vertical velocity from eq. {91 . 

The validity of the dispersion relation given by equation [9] is 
demonstrated by Figs. lA and IB, which show the vertical 
wavenumber as predicted by eq. [9] , denoted by open circles, in 
comparison with laboratory measurement (solid line) and the 
isotropic approximation given by eq. (10] (dashed line) in two 
perpendicular direction of the k space (k^. \) . 

Another form of q(z) can be defined by a transformation into 
polar co-ordinates ,i.e., writing q(z) as function of dip angle 
and azimuthal angle <|>. This can be done by changing the 
wavenumbers k^ and k, by the dip and azimuth angles: 



CO 



sin 9 cos (]). 



til] 





[12] 



k = — sin 9 sin ^ 



V 



With these new variables, the vertical wavenumber (equation [9]) 
can be written as a Fourier series in 4): 
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2 

ti3i q'(z) = % [Coce) + ;,(e) cos (2(t)) + C2<e) <=°^ (4<t»] - 



where 



10 



r (9) = 1 (v^ + ,_^) sin^ 9 - 

Sq^^' 2Y x.NMO y,NMO 

V 

i (2p^ + + 3Pg) sin* 6 
8 

[14) ^^(e) = - ^ <<NMO - <««d) Si"' ® - 

V 

i (Ps - Pg) sin" e 
;,(e) = i (2p, - Pj - p,)sin* e 

Equations [14] show that the orthorhoiribic vertical wavenximber 
can be described by three functions which are independent of 
azimuth (|>. The first part Co(9) represents an equivalent TI 
medixim. 

Using equations [13,14] it is possible to reconstruct the 
orthorhombic anisotropy by using three different three different 
azimuths. This is exploited in the scanning procedure for three 
specific values of azimuth ij> = 0, Ji/4, and JC/2: 



15 
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if<J> = 7t/2, with 

[18] k = ^1 + ^2 

5 These three expressions of q at three particular values of <J> 
conpletely describe the dispersion relationship, through the 
vertical wavenumber, in orthorhombic media. Furthermore, each of 
the three expressions describes a TI medium. Therefore, the 
orthorhombic velocity model can be parameterized by three TI 

10 media: For <j>=0, the vertical velocity, the normal moveout 

velocity and the anelliptic parameter of the corresponding TI 
medium are V^, ^x.nmo P5/ respectively. For the 

corresponding TI medium has the same vertical velocity but a 
different normal moveout, ^ and a different anelliptic 

15 parameter p6. The anellipticity of the TI medium corresponding to 
<j) = TC/4 contains the anelliptic paraoneter p4 needed to completely 
reconstruct the orthorhombic velocity model. 

In the following the above derived relations to reconstruct a 
20 background velocity model, including a method for selecting or 
sorting seismic data for such reconstruction, is described. 

As the anisotropic background velocity model is generally 
described by several elastic coefficients or anisot-ropy 
25 parcimeters, it is important to adopt an efficient scanning 

procedure, which effectively constrains these coefficients. In 
this example common azimuthal sections are selected from 
recorded three-dimensional (3D) seismic data. 



30 The procedure is then repeated for different common azimuthal 
sections, each new common azimuthal section leads to a new 



- 14 - 



azimuthally isotropic velocity model if the medium is 
azimuthally anisotropic. The number of common azimuthal 
sections, and therefore the number of azimuthally isotropic 
velocity models / needed to reconstruct an azimuthally 
5 anisotropic velocity model is dependent on the type of 

symmetries. For example, only three common azimuthal sections 
are needed for an orthorhombic medium. 

In the following a way of constructing a common azimuthal 
10 section from given seismic data is described: 

Considering a typical seismic reflection survey, (x^, y,, 2^=0) 
denotes a generic receiver position and (x^, y,, z^=0) denotes a 
generic "shot" position. A typical seismic reflection data sec 
15 can be represented by D(x,, y^, t; x,, y., ) . The time t is reset 
to zero at each shot . 

Introducing the midpoint co-ordinates 



X = 



(X + X ) / 2 



20 [19] 



y = 



and half-offset 



co-ordinates 



[20] 



h = 



h' = 



(x^ - x^) / 2 



25 



the data D(x,. y,, t; x., y., ) can equivalently be represented by 
D(x, y, t; h, h') or its Fourier transformed version D(k,, k^, 0); 
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K' K'^' The variables k^, k^, k,. and k/ are midpoint and half- 
offset wavenumbers which correspond to x, h^ and h' , 
respectively. 

5 To facilitate the scanning procedure for the reconstruction of 
azimuthally anisotropic background, thus, constant azimuthal 
sections are chosen instead of the full 3D seismic data. This is 
equivalent to taking one of the half offset co-ordinates, for 
exan^le h' , to be constant. A common azimuthal section can be 

10 described as a series of 2D multi-offset profiles all oriented 
according to a particular azimuth. For further details 
concerning the general aspects of this method, reference is made 
to L. Ikelle, Geophys . J . Int . , No. 123(1995), 507-528, and A. 
Canning and G.H.F. Gardner, 62nd Ann. Mtg. (1992), Soc. Expl . 

15 Geophys., Expanded Abstracts, 955-957, which describe how a 3D 
marine seismic accjuisition can be transformed into a series of 
2D profiles, from which the common azimuthal sections can be 
extracted. Savings with respect to CPU time and memory are 
gained from the use of constant azimuthal sections. 

20 

In practice, 3D multi-offset data are collected in multi- 
channels but not as a series of 2D multi-offset profiles. The 
sorting or organization of 3D multi-channel data into 2D 
profiles determines how a common azimuthal section must be 
2 5 defined. In this example it is simply assumed that common 

azimuthal sections are formed by brute stack over the offset in 
the direction perpendicular to the azimuth. For example, D(k,, 
k , CO; k^, k^'-O) • represents the zero degree azimuthal section. 

30 Alternative organizations from 3D multi-channel data to a series 
of 2D profiles are known. However the described method of 
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extracting a common azimuthal section from 3D seismic data 
simplifies the subsequent inversion formulae. 

Having selected the appropriate data, it is described in the 
5 following how to determine the optimal value for the anisotropy 
parameters and velocities as defined hereinabove. 

For the following steps it is assumed that a solution for the 
inverse problem, i.e., to determine the field description from 
10 the measured seismic data for an isotropic medium exists. Such a 
method is also described in L. Ikelle, Geophys . J . Int . , No. 
123(1995), 507-528, make use of the solution of the linearized 
inversion as known from L. Ikelle et al . Geophysics, Vol. 51, 
No. 6 (June 1986), 1266-1276, equation (14). 

15 

It is important to notice that alternative description for the 
Earth medium are known, leading to different solution of the 
inversion problem. For example a solution of the inverse problem 
in the space-time domain is given by A. Tarantola, Geophysical 
20 Prospecting 32(1984), 998-1015, summarized in equations (36). 

As all of these known solutions of the inverse problem are 
complex and use different notations, a reproduction of the 
relevant formulas in this description is omitted. 

25 

However, a common characteristic of the known solutions is that 
use is made of equations which characterize the movement of the 
wavefront of the elastic wave in space and time. In the f-k 
domain these equations are called dispersion relationships, in 
30 the space-time domain they are known as ^moveout-equat ions' , It 
is an important aspect of the invention to parameterize those 
relations, such that at least one dimensionless parameter 
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appears which characterizes the anisotropy of the medixim. 
Examples of such parameters are ri and p.. 

As a result of the parameterization of the dispersion 
5 relationship, the background medium is described by several 

parameters in the anisotropic case. Therefore, for setting up a 
scanning procedure in the anisotropic case, it is iiT?)ortant to 
scan over parameters that are as independent as possible, in 
terms of their relationship to travel times of seismic data, or 
10 over parameters that can be reconstructed hierarchically based 
on the significance of their contribution to travel times in 
seismic data. 

For reconstructing an azimuthally isotropic background velocity 

15 model (i.e., TI model), in the present example a two scanning 
step procedure is utilized. As the isotropic component is the 
dominant part for most anisotropic rock formations in the Earth, 
the first scan optimizes the normal moveout velocities, V^, 
which characterize the travel times variations with offset in 

20 isotropic media. This scan is carried out over isotropic 
background velocity models using linearized inversion in 
equation [1] with the dispersion relation in equation [4] . Once 
is determined, a second scan over the anelliptic parameter T| 
with the TI dispersion relationship of eq. [7] is carried out to 

25 adjust for anellipticity . Thus, the reconstruction of a TI 
velocity model requires two scans instead of one as in the 
isotropic case. In fact, an isotropic two-parameter linearized 
inversion, where the inversion of one parameter is essentially 
based on small offset data and the other parameter is 

30 essentially based on large offset data, yields, two different 
velocity models for the part of the data corresponding to 
anisotropic rock formations. One possible choice is the P-wave 
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impedance and the shear modulus with the appropriate weighting 
functions as chosen in this example (see above) . The normal 
moveout velocity in this case is extracted from outputs of P- 
wave impedance . 

5 

The optimum value for a given parameter can be determined from 
the scanning process performed in the above-mentioned steps 
using the various known criteria, such as Focusing criteria. 
Coherence criteria described by Yilmaz (see above) , which 

10 applies to Gaussian distributed data, or Bicoherence criteria, 
as known for exair^le from C.L. Nikias and J.M. Mendel, "Signal 
Processing with Higher-Order Spectra", United Signals & Systems, 
Inc., Massachusetts, 1990. Though the latter method gives 
probably the best results, in the present example the simple 

15 amplitude criteria is used, which selects the data according to 
their amplitude. 

Notice that the knowledge of vertical velocity is not required 
in this scanning procedure because the variations of vertical 
20 velocity do not affect the travel time 'variations with offset or 
azimuths; in other words, it acts as a scaling factor in time 
processing. However, the vertical velocity is crucial for depth 
conversion. Fortunately, it is often known from vertical seismic 
profiles (VSP) or check-shot data. 

25 

Having described above the process to optimize the parameter and 
velocity for a transversely isotropic medium, in the following 
the process of determining the optimal values for the parameters 
in case of an azimuthally anisotropic medium is described - 

30 

It has been shown that the dispersion relationship of 
azimuthally anisotropic media can be decomposed into TI 
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dispersion relationships (cf. equations [15]- [17] for the 
orthorhombic case. Therefore it is possible to reconstruct the 
azimuthally anisotropic velocity model as a group of TI velocity 
models . 

5 

Hence, an advantageous method for deriving the velocity model is 
to treat the problem as independent TI models. For the 
orthorhombic case, three TI model can be generated 
advantageously by selecting one azimuthal section of the data, 

10 and three different ranges of the azimuthal angle ^. For each 

range of (j) the two step scanning process yields one optimum value 
for a moveout velocity and one for an anisotropic parameter, 
which suffices to reconstruct the orthorhombic medium. The 
common azimuthal section is defined by the receiver positions, 

15 the azimuthal angle ij> describes the direction of the wavefront 
propaga t ion . 

Alternatively one fixed range of the azimuthal angle can be 
used, said range may extend to all possible values of <t>, together 
2 0 with three different sections of data with common azimuth. 

As another alternative, a combination of the two methods 
described before can be used to determine the optimal value of 
the moveout velocities and the anisotropic parameters of the 
25 medium. This combination is used in the following example, where 
the scanning process is repeated over different pairs of common 
azimuthal section and azimuthal angle. The steps are: 

(i) Selecting a zero degree azimuthal section and limit the 
30 linearized inversion to near zero degree azimuthal angle. 

Applying the two scanning step procedure to reconstruct a TI 
velocity model which is described in equation [15]. The 
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orthorhombic parameters thereby determined are the normal 
moveout velocity ^ and the anelliptic parameter . 

(ii) Selecting a 90 degree azimuthal section and limit the 

5 linearized inversion to near 90 degree azimuthal angle. Applying 
the two scanning step procedure to reconstruct a TI velocity 
model which is described in equation [17]. The orthorhombic 
parameters recovered here are the normal moveout velocity ^ 
and the anelliptic parameter p^. 

10 

(iii) Selecting a 45 degree azimuthal section and limit the 
linearized inversion to near 45 degree azimuthal angle, applying 
two scanning step procedure to reconstruct a TI velocity model 
which is described in equation [16], There is only one new 

15 orthorhombic parameter to be reconstructed here, the anelliptic 
parameter p^, because the normal moveout velocity in this common 
azimuthal section is the arithmetic average of the elliptic 
normal moveout velocities and V^^. 

20 Alternatively, zero degree azimuthal section or 90 degree 

azimuthal section without limitation on azimuthal angles can be 
used to adjust for the anelliptic parameter p^. The full 
orthorhombic dispersion [9] must be used in the linearized 
inversion. This alternative removes the need for a third common 

25 azimuthal ^ , V^^^ , p, and p, through steps (i) and (ii). 

In the following the above steps are carried out using a 
a geological model shown in Figure 3- The medium is composed of 
a homogeneous orthorhombic anisotropic material overlying an 
30 isotropic material. The interface between the two materials is . 
invariant with respect to the y-axis. The elastic coefficients 
of this orthorhombic material are (using the density normalized 
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stiffness tensor a(ijpq) in its compact notation A{IJ)) 
A{11}=12.79, A{22)=11.32, A{33}=8.56. A{12)=5.47, A{13)=5.14, 
A{23}=4.88, A{44)=2.30, A{55}=2.58, A{66}=2.76, V,=2.93, 
V,.^=3.47, V,._=3.22. p,=-0 . 27 , p,=0 . 0877 , p.=0 . 107 and by Figs . lA 
5 and IB, where use was made of the same medium. The A{IJ} are 
expressed in (km/s)' and the velocities in km/s; the anisotropy 
parameters are dimensionless . 

Two common azimuthal sections, zero degree azimuthal section and 
10 90 degree azimuthal section, were coir«>uted. The zero degree 
azimuthal section was formed as D(k., k,, a>; k,. k,'=0) and 90 
degree azimuthal section is formed as D(k., k,, <0; k,=0, k^'). For 
each common azimuthal section, 198 gathers corresponding to 
midpoints spaced every 12 m, with offset coverage between 100 m 
15 to 5000 m were computed using the Born forward modeling 

described by L. Ikelle, Geophys . J . Int . . No. 123(1995), 507-528 
with the exact vertical wavenumbers computed using the 
eigenvalue- eigenvector system. 
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AS stated above (step (i)) the process of background 
reconstruction is started with the zero degree azimuthal 
section. The inversion is limited to zero degree azimuthal 
angle. V,.^ is determined firstly by a scan over isotropic 
velocity models and then the anelliptic parameter p, is optimized 
by a second scan. The scan over velocities is carried out for 
velocities between 2900 m/s and 3700 m/s with an increment of 
100 m/s. The inverted parameter is the P-wave impedance. 
Applying the ait«)litude criteria to determine, an optimum 
velocity of 3500 m/s is established. The second scan is carried 
out over different values of anellipticity with V..^ = 3500 m/s, 
The anellipticity is varied from -0.1 to 0.3 with an increment 
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of 0.05. The best value for is between 0.05 and 0.1. The mean 
0.075 again is in good agreement with the actual value (0.0877). 



Then selecting the 90 degree azimuthal data. The inversion is 
5 limited to 90 degree azimuthal angle. Firstly, Vy.Nf« is 
determined by a scan over isotropic velocity models and 
subsequently the anelliptic paraimeter by a second scan. As 
result, the velocity could be determined as 3100 m/s or 32 00 
m/s, slightly below 

10 the actual value of 3220 m/s. This small discrepancy is due to 
missing small offsets. In fact, by working with a 90 degree 
azimuthal data restricted only to 90 degree azimuthal angle, and 
invariant with respect to the y-axis, the medium is inverted 
only at k^=ky=0, which equivalent to a one dimensional model. 

15 The missing near offsets have greater effect on the estimation 
of the normal moveout velocity when the medium is one 
dimensional . 

The second scan is carried out over different values of 
20 anellipticity with ^ = 3200 m/s. The anellipticity varies 

from -0.1 to 0.3 with an increment of 0.05. The best value of p^ 
is 0.10 which is the closest to the actual value for increments 
of 0.05. 

25 To completely reconstruct the orthorhombic background velocity, 
the anelliptic parameter p^ is determined. Using the 90 degree 
azimuthal data D(k^, k^, O); k^=0, k,^') without restriction on 
azimuthal angles and using the full orthorhombic dispersion, a 
final scan is carried out over different values of the 

30 anelliptic parameter p^. The values are taken between -0.2 and 0..2 
with an increment of 0.05 while the moveout velocities are fixed 
to 3500 m/s and 3200 m/s, and the anelliptic parameters p^ and p^ 
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are fixed to 0.7 5 and 0.10, respectively. The optimal value is 
determined as being -0.05, again very close to the actual value 
-0.027. 

5 Significant improvement in the imaging of dipping reflectors 
after correcting for the anellipticity is observed. 
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CLAIMS 

A 3-D seismic surveying method for surveying a subsurface 
earth formation comprising the steps 

- defining at least one moveout velocity and at least one 
anisotropy parameter; 

- inverting recorded seismic data using a plurality of 
different values for said velocity and said parameter in 

order to generate an optimized representation of said earth 
formation. 

The method of claim 1, coir^rising the step of determining a 
relation characterizing the movement of the wavefront of an 
elastic wave in time and space, said relation including the 
anisotropic parameter, and using said relation in the data 
inverting step. 

The method of claim 1, wherein an optimal representation is 
determined by using a criterion chosen from a group 
consisting of amplitude criteria, focusing criteria, 
coherence criteria, and bicoherence criteria. 

The method of claim 1, wherein the recorded data are prestack 
data. 

The method of claim 1, comprising the step of determining the 
representation of the earth formation by decomposing said 
representation into essentially independent representations 
of transversely isotropic media, each of said media having 
one velocity and one anisotropic parameter. 



6. 
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The method of claim 1, further including the step of 
arranging the seismic data into common azimuthal section. 



7. The method of claim 1, wherein the inversion is performed 
5 over a limited range of the azimuthal angle. 

8. The method of claim 1, wherein the number of anisotropy 
parauneters is less than 6. 

10 9. The method of claim 1, wherein the at least one anisotropy 
parameter is dimensionless . 

10. The method of claim 1 including the step of acquiring seismic 
data using means for generating waves travelling through 

15 subterranean formations and means to record said waves . 

11. The method of claim 1 including any combination of claims 2 
to 10. 

20 
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